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Abstract. We study the dynamic structure factor of a one-dimensional Bose gas 
confined in an optical lattice and modeled by the Bose-Hubbard Hamiltonian, using 
a variety of numerical and analytical approaches. The dynamic structure factor, 
experimentally measurable by Bragg spectroscopy, is studied in three relevant cases: in 
the clean regime, featuring either a superfluid or a Mott phase; and in the presence of 
two types of (quasi-)disordcred external potentials: a quasi-periodic potential obtained 
from a bichromatic superlattice and a box random disorder - both featuring a Bose 
glass phase. In the clean case, we show the emergence of a gapped doublon mode 
(corresponding to a repulsively bound state) for incommensurate filling, well separated 
from the low-energy acoustic mode. In the disordered case, we show that the dynamic 
structure factor provides a direct insight into the spatial structure of the excitations, 
unveiling their localized nature, which represents a fundamental signature of the 
Bose glass phase. Furthermore, it provides a clear fingerprint of the very nature of 
the localization mechanism which differs for the two kinds of disorder potentials we 
consider. In special cases, the dynamic structure factor may provide an estimate of the 
position of the localization transition from superfluid to Bose glass, in a complementary 
manner to the information deduced from the momentum distribution. 
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1. Introduction 

In interacting Bose fluids, the interplay between the effect of disorder and that of strong 
interactions displays a rich showcase of different behaviours. Ultracold atoms offer 
remarkably the possibility of exploring such an interplay, since in recent experiments 
both disorder and interactions are natural tuning knobs pQ. In these systems, disorder 
is realized either by the application of a speckle potential [21 El El El EE [7], orD y a 
bichromatic optical lattice composed of two incommensurate standing waves [HI El Ell 
CHI [12j H31 E2 US], giving rise to a quasi-periodic (QP) potential [16]. 

In a Bose fluid at zero temperature, Bose-Einstein condensation (or quasi- 
condensation in one dimension) occurs generically in the absence of disorder and in 
the weakly interacting limit. We consider two ways of destabilizing this phase on a 
lattice: either via Mott localization due to strong repulsion at commensurate filling; 
or via Anderson localization due to strong disorder at any filling [T71 IT8l [19]. The 
resulting phases - a Mott insulator for strong repulsion, and a Bose glass for strong 
disorder (and finite repulsion) - are very different forms of non-condensed Bose fluids. 
Yet, from the point of view of coherence properties, experimentally probed by time-of- 
flight measurements in cold-atom setups, the two phases look similar. In both cases 
phase correlations decay exponentially, giving rise to a broad peak in the momentum 
distribution. Hence additional information is required for a direct observation of the 
Bose-glass phase, which has been long sought in the context of Bose fluids. 

On the other hand, the Bose glass and the Mott insulator are fundamentally 
distinguished by the nature of their excitation spectrum. The lowest energy excitation 
in a Mott insulator is a particle-hole excitation with a gap imposed by the energy cost of 
the multiple occupation of a lattice site. In the limit of very strong repulsion, suppressing 
density fluctuations, such an excitation can be seen as a free particle/free hole moving 
on a static background of particles at integer filling. On the other hand, the lowest 
excitations in a Bose glass are gapless, and associated with phonon-like modes localized 
in rare, locally homogeneous regions of the sample. Therefore not only the density of 
states, but the spatial structure of the excitations provides a distinct fingerprint of the 
Bose glass phase with respect to the Mott insulating one. 

In this respect, the dynamic structure factor, probed by Bragg spectroscopy in cold- 
atom experiments [20], offers the possibility of characterizing both the spectral density 
of the excitations as well as their localized/delocalized nature. 

In the weakly-interacting regime, the dynamic structure factor can be estimated 
using the Bogoliubov approach [21], while at arbitrary interactions in the ID uniform 
system can be obtained from the Bethe-Ansatz solution of the integrable Lieb-Liniger 
model [22] [231 121], the long- wavelength behaviour in this regime is also captured by 
Luttinger liquid theory [25j. In the presence of a lattice, the dynamic structure factor has 
been the subject of several analytical (23 EH EH] and numerical [2S EDI ED E21 E31 M, ES] 
studies, as well as of recent experiments [2S1 EH EH1 ESI HD]- In the presence of a 
disordered potential, only a few studies have addressed the dynamic response functions 
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of lattice bosons such as the response for lattice modulation spectroscopy jS] and the 
single-particle spectral function |42], but we are not aware of previous studies of the 
dynamic structure factor. In particular the response to lattice modulation spectroscopy 
is sensitive to energy only, and it lacks the momentum resolution which, as anticipated 
in |43] from the study of the spatial Fourier spectrum of single-particle excitations, is 
essential to unveil the localized nature of excitations. 

In this paper, we aim at an extensive investigation of the dynamic structure 
factor in the case of one-dimensional Bose gases in an optical lattices, and for widely 
different regimes, encompassing the weakly interacting limit, the infinitely repulsive 
case (Tonks limit), and the regime of intermediate interactions. The two limits of 
weakly interacting and infinitely repulsive particles lend themselves to very convenient 
theoretical approaches (via Bogolyubov theory and fermionization, respectively), while 
the intermediate regime, which is the most challenging, can be investigated via exact 
diagonalization. We particularly focus our attention on the case of a quasi-periodic 
potential, featuring a localization transition for a finite potential depth even in the 
non-interacting limit - and we underline the analogies and fundamental differences with 
respect to a truly random potential. 

Our main result is that the dynamic structure factor serves as a very effective 
diagnostic tool of the localized phases, and, in selected cases, it might provide a 
quantitative method to estimate the quantum phase transition from superfluid to 
Bose glass, based upon the localization of the elementary excitations, and therefore 
complementary to the analysis of the coherence properties. In particular, the dynamic 
structure factor provides clear signatures of the underlying localization mechanism at 
play, in all interacting regimes, and it allows to characterize the Bose-glass phase far 
beyond its thermodynamic definition as a compressible and insulating phase. 

The paper is organized as follows. Section [2] introduces the model under 
investigation, i.e. the one-dimensional Bose-Hubbard model in a (quasi-) disordered 
potential, and the dynamic structure factor. Section[3]recalls some results on the limiting 
cases (weak interaction and infinite interaction) of the clean system, and it discusses 
exact diagonalization results interpolating between these limiting cases. Section [4] 
studies in details the disordered models: Section 4T describes the dynamic structure 
factor in the case of weakly interacting bosons, treated via a Bogolyubov approach, 
while Section 14.21 focuses on the exact solution in the case of hardcore bosons; Section 



4.3 bridges the two above regimes, contrasting the dynamic structure factor across the 



localization transition with the same quantity across the Mott transition in the absence 
of disorder. Section [5] contains a discussion on the relevance to experiments while Section 
[6] is dedicated to conclusions. 
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2. Bose-Hubbard Hamiltonians in the presence of disorder and methods 

2.1. Models 

We describe one-dimensional bosons in a deep lattice potential and subject to an external 
potential using the Bose-Hubbard Hamiltonian, 

n = - J J2fiU b 3 + h - c -] + |l>;(^ - !) + I>i^- • (!) 

Here bj is the operator creating a boson at site j, nj = b^bj is the local density, J is the 
hopping amplitude and U is the onsite repulsion. We consider lattices with L sites and 
N particles, i.e. a filling factor n = N/L. We denote the lattice spacing by a in the 
figures - when not appearing explicitly in the equations, it is understood that a = 1. 
The Wj are site-dependent energies which account for both the disorder potential and 
the possible presence of a harmonic trapping potential. For what concerns the disorder 
distributions, we will focus on two different forms of (quasi-) disorder: 

(i) A quasi-periodic (QP) potential obtained via a bichromatic optical lattice [8], with 
the form 

^ = ^^ 1 + cos ( 2r ^ + 2< ^ ( 2 ) 

j 

with r an irrational number and a random phase-shift on which averaging can be 
performed. We choose the experimentally relevant value r = 830/1076 [9]; when 
considering periodic boundary conditions on a lattice of size L, we take r as the 
best rational approximant in the form M/L (where M is a positive integer), so that 
the potential describes a single period over the entire lattice. 

(ii) A random-box (RB) disorder, for which Wj is a random variable uniformly 
distributed over the interval [0, V]. 

In both cases, V/J gives the relative strength of the (quasi-) disorder potential, 
and it is chosen in such a way that, in the atomic limit, the (quasi-) disorder closes 
the Mott gap for V — U. The main difference among the two types of disorder is 
that, in the absence of interactions, the QP potential leads to localization of the single- 
particle wavefunctions above a critical disorder threshold V = 4J [16] while for the RB 
disorder localization occurs at an infinitesimal disorder strength. The RB phase diagram 
was studied numerically in Refs. jlH H5j HHJ EH HE]. In the case of the QP potential, 
the localization mechanism and the features of the single-particle wavefunctions and 
spectrum was extensively studied [TB], HH1 EDI ETJ E21 E31 El]- In particular, a simple 
interpretation can be given in terms of successive band-folding processes (see e.g. 
Ref. [HU ESI S3]) or at a semi-classical level [SD ES]- Bosonization studies EH] 
also show that the quasi-periodic potential is different from the pure disorder one, in 
particular, it does not share the same expected universal Luttinger parameter value at 
the localization transition, which was checked in Ref. [13]. The phase diagram of the 
bichromatic system was investigated numerically in Refs. [591 EDJ H3J |6T] . 
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2.2. Dynamic structure factor 

The dynamic structure factor S(k, oS) is given by the space-time Fourier transform of 
the dynamic density-density correlation function, which for lattice bosons is given by 
(8rij(t)8ne(0)) with Snj = rij — (rij). It yields the linear response of the fluid to a density 
perturbation transferring a momentum hk and energy Hu to the system. Its Lehmann 
representation reads 

S(k,cu) = J2\^\n k \0)\ 2 5(u - u a ) , (3) 

where a labels the eigenstates of the Hamiltonian (0 being the ground-state), u a = 
E a — E Q are the excitation energies and 



<"'"■".' ( 4 ) 



is the density operator at momentum k, with rj = a(j — L/2). On finite lattice systems 
with periodic boundary conditions we use k = 2irm/L with integer m values. We will 
also consider the momentum integrated spectral function S{uj) = J2k S(k,u). 

The dynamic structure factor can be experimentally probed by Bragg spectroscopy, 
involving a two-photon transition. Using the fluctuation-dissipation theorem, it can be 
extracted from the measured energy gain of the system per unit time dE/dt, according 

to [621 EH] 

-—- oc u S(k, u) (5) 

where u is the frequency difference between the two photons involved in the transition 
and k is the wave-vector difference. Alternatively, the dynamic structure factor can be 
extracted from the rate of momentum transfer dPj dt [63] 

dP 

-j— oc k S(k, us) . (6) 

Both definitions have been exploited to extract S(k,u) in recent cold-atom experiments 
[HU EHl [65] [361 EHl HQ] . It is worth mentioning that the /-sum rule allows to compare 
theory with experiment without adjustable parameters [27] . 

Three different theoretical approaches to compute S(k, oj) are used throughout this 
work and we describe them below. 



2.3. Methods 

2.3.1. Exact diagonalization We computed S{k,u) using the Lanczos algorithm to 
represent the low-lying excited states. The method is exact but limited to small sizes. 
In the calculation, the maximum number of onsite bosons is fixed to 6, and 200 iterations 
are performed to compute the spectral weights. The delta functions in energy of the 
discretized excited states are convolved with lorentzians of width 0.2 J for S(k,u) and 
0.3 J for S{oj). Averaging is performed over 34 samples with a uniform distribution 
for the QP potential while 100 samples are used for the RB potential. 
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2.3.2. Bogolyubov theory For weakly interacting bosons having a finite condensate 
fraction in the ground state, it is well established that the ground-state properties 
as well as the low-energy excitation spectrum are well described within Bogolyubov 
theory. The traditional formulation of Bogolyubov theory requires the existence of a true 
condensate, and it amounts to neglecting terms which are not quadratic in the operators 
involving particles out of the condensate. Yet Ref. [66J has shown that an analogous 
approach to Bogolyubov theory can be applied to one- dimensional systems featuring 
only quasi-condensation. Such an approach is based on a polar decomposition of the 
Bose operators in terms of density and phase, bj = e*^*/nj, and on the fundamental 
assumption of weak quantum fluctuations of the density Srij = nj — pj around the mean 
Pj, as well as weak quantum fluctuations of the phase difference between neighboring 
sites 8j — Oj+i- The Hamiltonian can then be expanded in powers of the density and 
phase-difference fluctuations around a reference state, corresponding to non-fluctuating 
(and vanishing) phase differences and a classical density profile pj which satisfies a lattice 
Gross-Pitaevskii (GP) equation 

- J (y/Pj+i + y/PJ 1 ^) + [UPj ~ ~ w i)} VPj = ° ■ ( 7 ) 

Here p is the chemical potential controlling the number of bosons in the system 
(see below) and pj can be identified with the density profile of the quasi-condensate. 
Diagonalizing the quadratic Hamiltonian in the fluctuations amounts to a Bogolyubov 
transformation of the density and phase operators to operators a s , a\ 

= X/ {&n S j a s + 6n*j a\) + (d N pj) V 

s 

h = E fad a ° + a D - Q ( 8 ) 

s 

where V and Q are canonically conjugated operators associated with the zero-energy 
mode, and 



$n s>j = ^/p~ (u a j + v s>j ) 9 Sij = ' ' . (9) 

The u s j, v s j amplitudes satisfy the Bogolyubov-de Gennes (BdG) equations 

(10) 




V.2 Wu 





where 



U 2 = - jE(|j)(j + l|+h.c.) + E(2t/p J +^-/i)|j)(j| 

3 3 

Ku= ./)(/ (11) 

3 

and \u s ) = J2j u sj\j), \ v s) = Yuj v s,j\j)- The properties of the non-Hermitian eigenvalue 



problem of Eq. (10) are well known [S7|; in particular the solutions of the BdG equations 
with non-zero energy e s satisfy the normalization condition 

(u s \u' s ) - (v s \v' s ) = 8 SS > . (12) 
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In the absence of disorder the u and v modes are plane waves, and the corresponding 
energies, labeled by momentum (s — >■ k), have the well-known form 



\/e k (e k + 2nU) , with e k = 4Jsin 2 (A;/2) . (13) 



Given that quantum fluctuations in the density are linear in the a s , a\ operators, within 
the quadratic approximation for quantum fluctuations the quantum corrections to the 
density profile vanish (which is consistent with the image of pj as the mean local density). 
Therefore the total particle number is given by iV = Y2j Pji an d the chemical potential 
p in Eq. ^ is fixed so as to impose the desired lattice filling at the level of the GP 
equation. 

The Bogolyubov theory for quasi-condensates of Ref. [66J has been applied to ID 
Bose gases in a disordered potential in a number of recent papers (681 EHl El [71]. In 
particular Refs. [68], [691 E2] provide an explicit expression for the one-body density 
matrix, g^(j, I) = (b\b{) in terms of the solution of the GP and BdG equations, reading 



9 



( 1 \ J J) = (b}b l ) = v fpJpiexp 



It. 



V s,j V s,l 



y/Pl 



(14) 



where u^j, vj-j are the coefficients of the vectors \u s ) , \v s ) orthogonalized with respect 
to the quasi-condensate 

\uj) = \u s ) - (ip \u s )\ipo) (15) 

(and similarly for \v^r)). Here \ipo) = J^j^PojlJ) * s the normalized quasi-condensate 
mode, with ^Jp] = y/Nipoj. 

The calculation of the dynamic structure factor of a weakly interacting Bose gas 
has been addressed within Bogolyubov theory in various references [731 E] • m the case 
of quasi-condensates its expression turns out to be analogous to that of conventional 
Bogolyubov theory [74J, namely 

S(k,u) = Y,\SPs(k)\ 2 5(u-e s ) (16) 

s 

where the form factors Sp s (k) read 

8p s (k) = e " lkr ' 3 VPj( u s,j + v.j) • (17) 

3 

Notice that here we restrict our attention to the case u > 0, so that the zero- mode 
contributions disappear from S(k,u). In the absence of an external potential, the 
dynamic structure factor is a <5-peak resonance at the Bogolyubov mode e k 



S(k, u) = — N5( y co - u k ) . (18) 

The Bogolyubov theory for quasi-condensates is quantitatively consistent as long as 
the quantum fluctuations of the density and relative phase remain weak. In particular, 
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the phase-density formulation requires the condition pj ^> 1 to be fulfilled in order for 
the phase operator to be well defined as a (quasi-) Hermitian operator [66J. Moreover, 
a large quasi-condensate density on each site is also necessary for the relative particle 
fluctuations to be small, since ((Srij) 2 ) / p? > 1/pj [66]. As we will see, this condition 
will strongly limit the range of validity of our results. 

From the practical point of view, we numerically solve the Gross-Pitaevskii equation 
via split-operator imaginary-time propagation, and the Bogolyubov-de Gennes equations 
by diagonalization of the non-Hermitian C matrix using the LAPACK libraries, as 
described in previous references [75]. We present results for lattices with L = 256 and 
L = 512 sites. Unless otherwise specified, the results for QP potentials are averaged 
over ~ 50 values of the spatial phase <ft of the potential. 

2.3.3. Hardcore-boson limit While Bogolyubov theory applies to weakly interacting 
bosons at large filling, we can also consider the opposite limit of infinitely repulsive 
bosons, U — > oo, and low filling n < 1. This limit corresponds to the ID Tonks- 
Girardeau gas of hardcore bosons (HCB), in which the forbidden double occupancy of 
the sites can be incorporated in a redefinition of the bosonic operators, hi — > b{ satisfying 
bosonic commutation relations offsite and anticommutation relations onsite, \bj, b^] = 
for j 7^ / and {bj,bj} = 1, {bj,bj} = {b^b^} = 0. The hardcore boson operators bj can 
be transformed to fermionic ones Cj via a Jordan- Wigner transformation [76] , mapping 
exactly the hardcore boson Hamiltonian to free fermions with chemical potential p 

h = -jJ2 (sVi + h-c) + - ^ c h ( 19 ) 

which lends itself to efficient exact diagonalization. In the following we will focus on 
systems with open boundaries or in a trap, and therefore we omit in the Hamiltonian the 
boundary terms arising from the non-local nature of the Jordan- Wigner transformation. 

Most notably, the hardcore boson density coincides with the fermionic one, b\bj = 
CjCj, so that the dynamic structure factor of the hardcore bosons corresponds to that of 
the free fermions, taking the simple expression 

S{k,u) = M k )\ 2 f{ep,T)[l - f(e a ,T)} 6(u - cu af} ) (20) 

a/3 

where u a/ 3 = e a — e@, e a are the eigenenergies of the single particle problem in the QP 
potential, f(e,T) = {exp[(e — p)/{k B T)] + l} -1 is the Fermi-Dirac occupation factor at 
a finite temperature T, and 

Pa p{k) = Y J ^r a , j ^- (21) 

3 

Hence we observe that, for hardcore bosons, the fc-dependence of the dynamic structure 
factor describes the power spectrum in momentum space of the overlap function ip* a j^pj 
between occupied and unoccupied single-particle states, connected by the energy transfer 
ftw. 
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3. Dynamic structure factor of the clean system: superfiuid and Mott 
phases 



The phase diagram of the clean Bose- Hubbard model displays two phases [19j : the 
superfiuid phase (SF) which occurs generically at incommensurate densities, and the 
Mott-insulator (MI) phase which occurs only at commensurate fillings and beyond 
a critical interaction strength U c (for a one-dimensional system with filling n = 1, 
U c ~ 3. 3 J). In the incommensurate case at densities n < 1, there are two main types of 
low-energy excitations contributing to the dynamic structure factor in the Bose-Hubbard 
model. The first type is represented by gapless acoustic modes related to the superfiuid 
regime and which have, in the long wavelength limit, the dispersion relation u(k) ~ ck, 
with c the sound velocity. The second type is represented by doublon excitations, namely 
repulsively bound states of two particles occupying the same site, occurring when the 
repulsion energy exceeds the bandwidth; these bound states appear at the two-body 
level, and survive in the many-body case [77J . Its energy creation cost is about 4 J + U 
in the strong coupling limit. In the Mott regime, the elementary excitations form a 
particle-hole continuum. This gapped excitations of the Mott phase have a typical 
energy cost of the order of U and the shape of the particle-hole continuum is known in 
the deep Mott limit (see e.g. Ref. [2TJ ) . This excitation is essential in understanding 
the dynamical properties of the system at large U, and it is the main excitation in the 
Mott phase where sound modes are absent. 

These elementary excitations are visible in the dynamic structure factor of the clean 
Bose-Hubbard model. We give in Fig. [ija) the evolution of the dynamic structure factor 
S(k,u) for increasing interaction U/J and for two typical densities: incommensurate 
(n = 0.5) and commensurate (n = 1). 

We start with the incommensurate case. At small U / J, we display the Bogolyubov 



result of Eq. (18). The overall behavior has a form of an arc on [0, 2ir] due to the 
periodicity in momentum space. For U = 2J, we see that the Bogolyubov result 
and the exact diagonalization result differ qualitatively. While it is expected that 
Bogolyubov theory fails to account for such a relatively strong interaction, we observe 
that the essential difference emerges around k ~ 7r where local physics is dominant. The 
spectrum is there split into two lines which we interpret as an hybridization between 
the acoustic modes and the gapped doublon state. Then, the peak emerging on top of 
the acoustic branch is attributed to the doublon state. This doublon mode is roughly 
centered around the energy 4 J + U, as we see for increasing U/J, and as expected for 
repulsively bound pairs [77] . When its energy increases, its spectral weight decreases, 
as one can see from S(co) plotted in Fig. W[b). In addition, increasing the interaction 
strength transfers spectral weight to the low-energy states at the k = 2k F wave vector, 
k = 2im {i.e. k = 7r in the figure) corresponding to back-scattering (see eg [78]). This 
is the very analogue of what is found in the Lieb-Liniger model [22] in the absence of 
the lattice. In the hardcore boson limit, the spectrum corresponds to the XX spin chain 
model which displays the famous Pearson-De Cloiseaux continuum (see e.g. [791 El]). 
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Figure 1. (a) Dynamic structure factor S(k,w) for the clean system. Upper panel: 
for an incommensurate filling n — 0.5 (L = 22) from weak to strong interactions. 
Bogolyubov and hardcore-boson (HCB) results are compared to exact diagonalization. 
Lower panel: for a commensurate filling n = 1 (L = 16) where the Mott phase emerges 
at large U/J. The perturbative (pert.) result is also shown. Color scales from different 
figures arc different, (b) Integrated dynamic structure factor S(lj), allowing to compare 
the spectral weights for different U / J for the two situations. 



In this limit, the doublon excitation is no longer part of the spectrum. Thus, exact 
diagonalization nicely interpolates between the two regimes, and it exhibits the evolution 
of the doublon excitation in the spectrum when approaching the HCB limit. 

Looking at the commensurate n = 1 case in Fig. [TJ one can see the transition 
from the SF to the MI phase when interactions are increased, although the opening 
of the gap in the dispersion relation is appreciable only at U sizably larger than U c 
because of finite-size effects (the gap opens exponentially slowly). One switches from 
the SF behavior at U = J, similar to the incommensurate case and in qualitative 
agreement with Bogolyubov theory, to a fully gapped excitation spectrum in the MI 
regime which corresponds to the particle-hole dispersion centered around u ~ U. For 
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U = 2,7, there is a single peak in the acoustic mode, with a short lifetime at k = n 
compatible with previous results for the same filling [33j EH [35] . Interestingly, in the 
intermediate interaction regime U ~ 6J, the system is gapped but the particle-hole 
dispersion has a significant weight at the lowest frequencies, reminiscent of the acoustic 
mode spectrum [2B]- According to perturbation theory [27J, for large enough U/J the 
spectrum is predicted to display a butterfly-like shape with a maximum weight around 
k = 7T. For U = 10 J, where the gap is sizeable, S(k,u) has a similar support as the 
perturbative prediction of Ref. [27], but the weight distribution is not yet symmetric 
around u = U. We observe that ED can quantitatively covers the evolution of S(k, u) 
from the weakly to the strongly interacting regime, and that Bragg spectroscopy can 
capture the opening of the gap. Yet, in realistic conditions with a trap, one would have 
to consider the effect of the inhomogeneity of the system [27J. 

4. Dynamic structure factor for the quasi-periodic system 

The impact of disorder on elementary excitations of interacting bosons can be expected 
to be qualitatively similar to that on single-particle states, leading in particular to 
localization of the spatial support of the excitation modes. The connection between 
single-particle and many-body physics is evident within the Bogolyubov and HCB 
approaches, since the spatial structure of the excitation modes comes from the solution 
of the single-particle Schrodinger's equation (for HCB) or the solution of BdG equations 
in the presence of a (quasi-) disordered potential. Yet the same connection is far less 
obvious in the full Bose-Hubbard model. In what follows we will first describe our 
results in the Bogolyubov and HCB regimes, and then show how exact diagonalization 
allows to interpolate between the above regimes. In the case of exact diagonalization we 
also compare the case of QP and RB potentials, and we show that, due to the different 
nature of the localization mechanism at play, these two potentials lead to very distinct 
features in the dynamic structure factor. 

4-1. Results from Bogolyubov theory 

We present here the results for the dynamic structure factor of one-dimensional weakly 
interacting bosons in a quasi-periodic lattice. We begin our discussion with the non- 
interacting limit, which serves as a useful reference for the results in the interacting 
case. 

4-1.1. U=0 In the case of an ideal gas, the dynamic structure factor takes the simple 
form 

S(k,uj) = J2\P a o(k)\ 2 5(u J -Lo a ) (22) 

Q 

where J2 a runs over the single-particle eigenstates, u a = e a — e is the excitation energy 
of the a state, and e a is the single-particle eigenenergy corresponding to the lattice 
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Figure 2. Dynamic structure factor S(k,u>) for free bosons in a quasi-periodic 
potential for increasing strength V/J of the disorder potential. The case V = 2 J 
is shown in log-scale to highlight the underlying excitations bands. 



eigenfunction ip a ,j] Pao{k) is the Fourier transform of the overlap function between the 



ground-state and the a-th excited state, defined in Eq. (21). 

If the ground-state is close to a k = plane wave, the p a o{k) form factor is essentially 
proportional to the Fourier transform of the excited state ip a . As a consequence S(k,u) 
gives the power spectrum in momentum space for the excited state at energy e a = u. 
Therefore, even in the absence of translational invariance (broken by the QP potential), 
the presence of a sharp ridge in S(k, u) gives an effective energy-momentum dispersion 
relation for the excited states of the system. 

Fig. [2] shows the single-particle dynamic structure factor for an increasing strength 
of the QP potential. We observe that the e a = dispersion relation of free particles 
in the lattice, characterized by a single cosine band, breaks up into sub-bands for a 
finite V. We shall first focus on the delocalized phase V < V c . The appearance of the 
sub-bands can be related with the fact that the QP potential introduces a quasi-periodic 
structure in the lattice, whose spatial period corresponds to the period of the beating 
between the underlying lattice and the incommensurate potential, /qp = (1 — r) -1 . 
Correspondingly features in momentum space appear at the edges of a (pseudo-)Brillouin 
zone with width /cqp = 27r(l — r), namely at /cqp/2, 7t±/cqp/2, etc.; these are indeed the 
(approximate) momentum locations at which the gaps between the sub-bands appear 
in S(k,u). Within each subband the excitations are delocalized with sharp momentum 
content, and they exhibit a cosine-like dispersion with the periodicity of the pseudo- 
Brillouin zone. Yet, due to the incommensurability, the pseudo-Brillouin zone cannot 
fill the Brillouin zone of the underlying lattice an integer number of times, and hence 
the sub-bands dispersion curves essentially fade away in S(k,u) after a few periods (in 
fact a closer inspection shows that they persist over the whole Brillouin zone, and they 
even wind around it giving rise to a very complex pattern, which nonetheless is only 
seen in logarithmic scale - see Fig. [2] for V = 2J). This fading dispersion relation 
can be understood within a perturbative picture for the QP potential: a particle with 
momentum k and energy e& (in the absence of the QP potential) is scattered by the 
QP potential and it can acquire a momentum p /cqp at p-th order in perturbation 
theory, but due to the incommensurability there is no finite order in perturbation theory 
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Figure 3. One-body density matrix for the ID weakly interacting lattice Bose gas in 
a QP potential of increasing strength. The parameters are r = 395/512, U = Q.1J, 
n = 10, and L = 512. 



which can connect the initial state to a resonant unperturbed state, and therefore the 
particle remains "localized" around its initial momentum k with fading components at 
k ± k ± 2/cqp, etc. This picture of localization in momentum space is valid beyond 
perturbation theory, and it relies on the exact duality of the non-interacting model under 
Fourier transformation [161 [50] . 

For V > V c , S(k, us) undergoes a radical change: the dispersive nature of the 
excitations within the subbands disappears, and S(k, uS) acquires features which are very 
broad in momentum space, while retaining a sharp nature in the frequency domain. This 
corresponds to the appearance of strongly localized modes, possessing a large uncertainty 
in momentum space. The large broadening of the structure of S(k, u) in momentum 
space is therefore the signature of localization, and it will reappear as a leitmotiv in the 
analysis of the results for the interacting system. 

4-1-2. U > 0. In the following we present our results for the weakly interacting case. 
We will mostly present results for U — 0.1J and a lattice filling n — 10 to satisfy the 
conditions of validity of Bogolyubov theory - although we have also investigated the 
interaction strengths U = 0.01J and 0.5 J, displaying similar features to U = 0.1J. As 
already discussed in Refs. [HOI ESI ESI EE] for the case of bosons in continuum space, 
Bogolyubov theory is capable of describing quasi-condensates, and specifically a power- 
law decaying one-body density matrix. This is also verified for a lattice system and 
in the presence of a QP potential, as shown in Fig. [3j In particular we observe that 
weak interactions make the quasi-condensate state robust to the QP potential, and 
they promote it to values of V well beyond the critical value V c = 4J for the non- 
interacting system. In the case of a ID gas in continuum space and subject to a speckle 
or quasi-periodic potential, Refs. [BH1 ESI EI] have shown that Bogolyubov theory allows 
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to quantitatively describe the localization transition in one dimension in the presence 
of interactions; such a transition is detected by the appearance of an exponential decay 
in the one-body density matrix. This has allowed the authors of Refs. [6U [691 E] to 
track the interaction-induced shift of the critical disorder strength. In the lattice system 
under investigation, on the other hand, we rather find that Bogolyubov theory fails to 
reproduce quantitatively this transition. Indeed, for all the interaction strengths we 
considered, we find that a quasi-condensate phase is observed over the whole range 
of applicability of the theory, namely for disorder strengths V which do not lead to 
excessive fragmentation of the density profile. Indeed, if the disorder strength is too 
large, there will appear sites in the lattice with -C 1, clearly violating the condition of 
weak density fluctuations. We do observe a change of the g^ 1 ' function from an algebraic 
to an exponential decay with increasing V, but this occurs at unrealistically large values 
of V, well beyond the value ~ V c + Un which naively represents the critical value for a 
QP potential screened by the interactions, and well beyond the range of validity of the 
theory. 

Even if Bogolyubov theory does not allow us to describe the localized Bose glass 
phase for the ground state of the system under investigation, it still reveals a dramatic 
evolution in the properties of the excitations, and a very peculiar nature of the persistent 
quasi-condensate phase protected by the interactions. The evolution of the dynamic 
structure factor for an increasing strength of the QP potential is shown in Fig. |4j When 
comparing it to the non-interacting case of Fig. [2j one clearly observes substantial 
analogies. In particular under the effect of the QP potential, Bogolyubov modes are 
still organized in sub-bands, which exhibit sharp dispersion relations in the (k, oj) plane 
for sufficiently weak V, while they lose completely their definition in momentum space 
when the modes undergo localization for a larger value of V. We observe that the 
modes at higher energy localize at a lower value of V, as we will further elaborate upon 
in the following. In particular the lowest sub-band, containing the gapless excitation 
modes above the ground state, preserves its dispersive nature for all the values of the QP 
potential considered, even if the bandwidth gradually decreases with V - this is exhibited 
in Fig. [4] for V = 16 J, where a low-energy zoom on S(k,co) is presented for the two 
strongest values of V shown in Fig. |4} In particular the effective dispersion relation of 
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Figure 5. Sound velocity as a function of the QP potential strength and for various 
values of the inter-particle repulsion. Other parameters as in Fig. [4] The dashed curves 
are parabolic fits c(V) — c(0) — jV 2 for the lowest V values. 



the lowest band preserves a linear behavior for k — >■ 0, characteristic of a delocalized 
sound mode. From the slope at k — > 0, we extract an effective sound velocity c, which 
is shown in Fig. [HJ We find that c decreases as c(V) ~ c(0) — ^(U, n) V 2 where 7 
is a constant; this is consistent with the perturbative results of Ref. [E21 E3], showing 
that the quadratic dependence on V is a generic property of Bogolyubov modes in the 
presence of an external scattering potential. 

The progressive localization of Bogolyubov modes from the higher to the lower 
energies upon increasing the QP potential can be quantitatively captured by inspecting 
the effective spatial support of the u and v lattice functions. Following the natural 



definition of norm for the it, v functions as in Eq. (12), and in analogy with the case 
of normalized wave functions, one can define a participation ratio (PR) for the u, v 
functions in the form 



PR,,,, (a; 



1 2 

1 1 



L Ej «,■ - '•>.,) 2 L Ej ("'-,, - d if 



(23) 



This quantity captures the fraction of the system size over which the it, v functions 
have a non-negligible value. Fig. [6]^a) shows the evolution of PR™ as a function of both 
QP potential strength and excitation energy: it is clear that for every finite value of V 
the high-energy Bogolyubov modes are more localized than the low-energy ones, and in 
particular for V ~ 6J the highest band of Bogolyubov modes undergoes localization, 
and the lower bands follow in the localization cascade at higher V, while the lowest 
band remains delocalized over the entire range of V values covered by the figure. 



Coming back to the dynamic structure factor, its expression, Eq. (16), probes 
the spatial structure of the overlap function Jp~j(u s j + v 8 J) - giving the local overlap 
between the quasi-condensate mode and the excitation mode - and not simply the spatial 
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structure of the u, v functions. Nonetheless, if the condensate mode is delocalized, then 
the overlap function has the same localization properties as the u, v functions. This can 
be directly inspected by plotting the participation ratio for the overlap function: 



PR-™o(w) 



s,3j 



L Pi ( u s,j + v, 



\2 



(24) 



As shown in Fig. |6|a), the behavior of this quantity (as a function of QP potential 
strength and mode energy) is qualitatively very similar to the participation ratio of the 



u, v functions of Eq. (23). As a consequence, the dynamic structure factor can capture 



the localization properties of the u, v functions given that in its expression, Eq. (16), 
the Bogolyubov modes at energy e s are weighted by the power spectrum \Sp s (k)\ 2 of the 
overlap function. In particular the power spectrum has opposite localization properties 
with respect to the overlap function, namely it is delocalized in k space when the overlap 
function is localized and vice versa. Therefore, it appears natural that the localization 
properties of the overlap function can be extracted from the dynamic structure factor 
by examining its inverse participation ratio (IPR) of in k space 

IPRsH = • (25) 

In particular one can easily show that S(oj) = LJ2jPj( u s,j+v s ,j) 2 , so that IPR5 
and PR u „o share the same denominator. Fig. ^b) shows IPRs as a function of V 
and u; a comparison with the (V, oS) dependence of PR U do shows striking similarities, 
demonstrating that the dynamic structure factor allows to measure directly the 
localization properties of the excitation modes in the Bogolyubov regime. 

4-2. Results for hardcore bosons 

4-2.1. Weak QP potential. For a weak QP potential the (pseudo-) dispersion relation of 



single particles is altered as discussed in Section 4.1.1 with the opening of gaps at /cqp/2, 
7r ± kqp/2, etc... In Fig. [7] we show the dynamic structure factor for hardcore bosons of 
variable density in a weak QP potential V = 0.25J. In one- dimensional free fermions, 
some of the dominant features in the structure factor are related to the transitions 
between the states at the bottom of the dispersion relation and states anywhere else in 
the energy spectrum - this is due to the singular contribution of the low-energy states, 
associated with the van-Hove singularity in their density of states. As a consequence we 
observe all single-particle gaps in the dynamic structure factor as long as Pauli principle 
allows the corresponding transitions, namely as long as the Fermi wave vector, kp = nn, 
is lower than the wave vector of the arrival state - otherwise the transition is forbidden by 
Pauli blocking. Indeed we see that the gap at k ~ /cqp/2 disappears when k F > kq P /2 
(corresponding here to n > 0.23), and the low- A; and low-o; structure factor is dominated 
by the linear mode with dispersion 2Jsin(kp) k. In what follows we will focus on this 
situation, and investigate the case of (local) filling n m 0.3. 
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Figure 6. (a) Participation ratios PR ut , and PR u „o ( see text for the definition) of 
variable V and uj. Here U — O.f, r = f98/256 and L — 256; the data are referred to 
one single realization of the phase of the QP potential. The energy axis is discretized 
into intervals of width 5lo = 0.02 J, and the participation ratios (PRs) shown in the 
figure are an average of the values of the PR for the various modes falling within the 
energy interval. The dark-blue regions - associated with a vanishing PR - correspond 
to the gaps between Bogolyubov bands, (b) Inverse participation ratio of the dynamic 
structure factor; same parameters as in the previous panel. 
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Figure 7. Dynamic structure factor for hardcore bosons in a QP potential with 
strength V = 0.25J. Here we consider an open chain with L = 200 and variable 
density, and a single realization of the QP potential. 




Figure 8. T = dynamic structure factor of ID hardcore bosons in a QP potential 
with variable strength, and in the presence of a confining parabolic potential. 



4-2.2. Strong QP potential and localization transition. In the presence of a QP 
potential with strength V > 4J all single-particle states localize, and therefore we 
expect a radical change in the k dependence of the form factors p a p{k) in Eq. 20, as 
the overlap functions ipaj^pj ev °l v e from extended to localized. This is indeed observed 
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in Fig. [8] where we consider the evolution of the dynamic structure factor across the 
localization transition at T = for a system of N = 60 hardcore bosons in a QP 
potential, and further confined by a weak harmonic potential Vjrj (to make contact with 
an experimentally realistic situation) with Vt = 10~ 3 J. We observe that the dispersive 
nature of the excitation modes is quickly lost as V approaches the critical value, and that 
the structure factor fragments into horizontal ridges, namely excitation modes which are 
well defined in energy but poorly defined in momentum space. Such features correspond 
predominantly to localized particle-hole excitations, in which the two states ip a and ipp 
connected by the transition are both localized in the same region of the system, giving 
a sizable overlap function. 



4-3. Exact diagonalization results: competition between Mott insulator and Bose glass, 
and comparison between QP and RB potentials 

4-3.1. Exact diagonalization results for QP potentials We now turn to the results of 
exact diagonalization in the QP case. The dynamic structure factor is given in Fig. [9] 
for four typical situations with increasing QP potential: incommensurate with weak 
and strong interactions (Fig. (9^a)), and commensurate for the same interaction regimes 
(Fig. [9^b)). In the incommensurate case, the effect of the bichromatic potential is 
rather weak at U = 2J and V = 2J. The u width increases and the doublon mode 
is hardly visible but one does not see the opening of gaps expected from Bogolyubov 
theory, possibly because of finite size effects. Increasing further the potential leads to 
localization of the excitations and a spectrum qualitatively similar to the one predicted 
using the Bogolyubov approach in Sec. |4.1.2 Notice that although the spectrum 



apparently looks gapped due to a large spectral weight for excitations at ui ~ 3J, it 
is in fact gapless but with small weights for low-energy excited states. The effect of the 
QP potential is more evident when starting from U = 10 J (close to the HCB limit) and 
increasing V. Subbands in the spectrum do appear (see the panel with U = 10 J, V = 2 J 
in Fig. [9]) while the doublon mode loses some of its dispersive features, consistently with 
the fact that it can be localized already at weak QP potential strength - it has a reduced 
effective hopping ~ J 2 /U - and its energy is lowered by the disorder. For a strong QP 
potential (V = U), the spectrum exhibits many subgaps and it strongly broadens in 



k, in a very similar manner to what seen in the HCB results of Sec, 4.2 The doublon 
mode is no longer resolved and the spectrum has large weights over a broad range of 
frequencies, while the integrated weight S(u>) shows a strong suppression. 

Turning to the commensurate case in the superfluid regime (U = 2J), here again 
a sufficiently large bichromatic potential is required to change the spectrum. Above 
the localization transition, the spectrum displays many subgaps typical of the band- 
folding localization mechanism (as seen in the panel for U = 2 J, V = 10 J in Fig. gb)) 
before reaching similar strongly localized spectrum as in the incommensurate density at 
very large V (see the V = 20 J panel). Starting from the MI phase at U = 10 J, and 
introducing a weak QP potential (V = 2 J in Fig. gb)), we observe that the spectral gap 
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Figure 9. Dynamic density structure factor for the bichromatic system, (a) for a 
typical incommensurate density n = 0.5 (L = 20). (b) for the commensurate density 
n = 1 (L = 14) where the Mott phase emerges at large U/J. 

is initially lowered, and mini-gaps appears in the particle-hole dispersion. For a stronger 
QP potential (V — U — 10 J in Fig. [9](b) ) , the gap closes and the system enters the 
strongly-correlated Bose glass phase. The spectrum exhibits both low-energy excitations 
with weights around k = - corresponding to phonon-like modes of locally superfluid 
regions - and excitations at relatively high energies - corresponding to short-wavelength 
localized excitations. Increasing further the QP potential to largely exceed the Mott 
gap (V = 20 J = 2U in Fig. [9^b)), the spectrum appears as composed of two parts: 
a low-energy part, associated with regions exhibiting locally incommensurate densities, 
and hence a similar behavior to that of the incommensurately filled lattice; and a higher- 
energy part with uj ~ U, associated with localized particle-hole excitations appearing in 
regions with local Mott behavior at commensurate filling. A similar separation emerges 
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with the RB distribution as we will see in the following section. As an intermediate 
conclusion, one can keep in mind that the typical signature of the localization due to 
the QP potential is best observed when the spectrum possesses subbands and is very- 
broad in k. The spectra for U = 2J, V = 10J, e.g., displays nicely this fingerprint. 



4-3.2. Comparison with random box disorder We now compare the effect of a QP 
potential seen in the previous section with the effect of true disorder, represented by the 
RB distribution. We show results for a RB potential having the exact same values of the 
strength V as those discussed for the QP potential. As already discussed in Ref. [43J, the 
Fourier transform of single-particle excitations for a RB potential differs significantly 
from those for a QP potential. In particular the gaps occurring for a QP potential are 
absent for a RB potential, and momentum broadening due to localization occurs at an 
infinitesimal strength of the RB potential, while a weak QP potential rather leads to 
excitations with narrow momentum features, and with a quasi-period imposed by the 
pseudo-Brillouin zone of the QP potential. This is due to the fact that a RB potential 
scatters Bloch waves at any wave vector, while a (weak) QP potential primarily affects 
Bloch waves with a wave vector k ~ kqp/2, 7t±/cqp/2, etc. On the contrary, a strong QP 
potential induces a stronger localization than the RB, leading to very broad momentum 
features. These differences between the RB case and the QP one will clearly persist in 
the Bogolyubov and HCB regimes. In what follows we check that similar differences 
are also present for interaction strengths interpolating between the weakly and strongly 
interacting regimes. 



The results for S(k,u) in the presence of a RB potential are plotted in Fig. 10 
The values for the interaction and disorder strengths are the same as in Fig[9j To avoid 
repetitions with respect to the discussion of the QP case, we simply highlight the main 
analogies and differences between the RB and QP case. For the incommensurate density 
n = 0.5 and weak interaction U = 2 J, we generally observe that an arc-like shape of the 
S(k,uj) support, typical of the Bogolyubov mode in the clean case, is preserved in the 
presence of disorder, but the energy and momentum structure is strongly broadened, 
due to localization of the modes and to the random distribution of energies of localized 
excitations induced by disorder. This is in sharp contrast with the subband formation 
seen in the QP potential. The doublon excitation undergoes a similar fate as for the 
QP potential, merging rapidly with the acoustic modes for a large enough RB potential. 
For a stronger repulsion U = 10 J the evolution of S(k, u>) with increasing disorder is 
comparable to the quasi-periodic case, certainly due to the fact that finite-size effects 
make the small differences hardly visible. Yet, in the large disorder limit, the spectrum 
is quite different from the QP, with the absence of subbands and a weight broadly 
distributed in frequency. In both cases, the spectral weight tends to decrease with V. 

Turning to the commensurate case n — 1 at weak interaction U — 2 J, we find that 
the spectra at weak disorder are comparable to the QP ones, the disorder leading just 
to broadening of the dispersion, while the spectra at strong RB disorder are remarkably 
different from the QP case, illustrating the different localization mechanism in the two 
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(a) incommensurate n = 0.5 with box disorder 




Figure 10. Dynamic density structure factor for the RB disordered system, (a) for a 
typical incommensurate density n = 0.5 (L = 20). (b) for the commensurate density 
n=l(L = 14). 

cases. In the case of strong interactions, U = 10 J, we observe that the RB potential leads 
to a closing of the Mott gap, similarly to the QP potential, but deep in the Bose-glass 
phase (V = 20 J), S(k, u>) exhibits a rather special structure with two coexisting features: 
a low-energy arc-shaped part, quantitatively consistent with the incommensurate filling 
case at the same disorder strength (compare the picture at U — 2J), and a high-energy 
part with the same structure as the particle-hole excitations of the Mott insulator at 
weak disorder (compare the case V = 0.4 J). As for the QP potential, this is a clear 
signature of the strongly-correlated Bose glass regime, with the coexistence of regions 
with locally incommensurate filling and gapless excitations, and regions which preserve 
a commensurate filling and a Mott-like behavior - as also seen (albeit less clearly) in 
the case of the QP potential for the same interaction and potential strength. 
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Figure 11. Dynamic structure factor of ID bosons in a QP potential with variable 
strength and in a confining parabolic potential with strength V t = 10~ 3 J. Here the 
QP potential is increased adiabatically at fixed entropy per particle S/N = 1 k^. 



An important conclusion of this section is that, the comparison between QP and 
RB results at strong disorder shows that Bragg spectroscopy could certainly probe the 
very nature of the Bose glass phase and unveil the localization mechanism at play. 



5. Experimental considerations: finite entropy and w-scan overlaps 

In the previous section we have seen that low-energy features of the spectral function 
at large U reproduce the HCB behavior. Therefore, one can exploit the exact solution 
available for the HCB case to make further contact with an experimentally realistic 
situation. Given that in experiments the loading of the optical lattices occurs in a 
(quasi)-adiabatic way, we study the evolution of the dynamic structure factor with an 
increasing height of the secondary lattice at fixed entropy per particle, taken to be 
S/N = Iks- The corresponding temperature which enters in the calculation of the 



dynamic structure factor in Eq. (20) is obtained by numerical inversion of the equation 



which links temperature and entropy for free fermions 

S(T) = -k B J2 [fa ^g f a + (1 - /„) l0g(l - fcx ) ] (26) 

a 

where f a = f(e a ,T). In the present case, since the QP potential reduces the density 
of states at low energy and fragments the energy spectrum into increasingly spaced 
minibands, entropy conservation implies adiabatic heating of the system as V increases. 
Yet, the comparison between Figs. [8] and 11 shows that the main features of the 



localization transition remain intact, indicating that they are accessible to current 
experiments. 

As a second aspect with direct experimental relevance, we propose an effective, 
global way to capture the dispersive or non-dispersive (namely fc-dependent or k- 
independent) nature of the dynamic structure factor. This amounts to considering 
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Figure 12. (a) Overlap function of w-scans in the dynamic structure factor for N = 60 
hardcore bosons at T = and in a QP potential plus a confining parabolic potential. 
Each panel shows the overlap at four different reference wave vectors fc, with an 
increasing separation Ak when going from the upper-left to the lower-right panel, 
(b) Same but for finite entropy, S/N = 1 kg- 



the overlap function of two w-scans in S(k, u) at wave vectors k and k + Ak: 

/ duj S{k, u)S(k + Ak, u) 

0{k,k + Ak) = 



1/2 

dco S 2 (k,oo) I dco S 2 (k + Ak,u) 



(27) 



This expression is normalized so that 0(k, k) = 1, and in general 0(k, k + Ak) is close 
to 1 if the features in S(k, u) at k and k + Ak have a large overlap; if on the other hand 
S(k, uj) has strongly dispersive features, O is in general <C 1 - it exactly vanishes in the 
extreme limit of a S{k,u) ~ 5{u — e^), displaying a 5-ridge associated with a dispersion 
relation having a finite group velocity around k, de^/dk ^ 0. 

In Figs, p^a-b) we show the overlap function at zero and finite entropy (1 
per particle) respectively, for four different values of k (going from the Brillouin zone 
center towards the edge) and for four values of the separation Ak between cu-scans, as 
a function of the QP potential strength V. For all scans we observe that the overlap 
increases for V/J — > 4™, and it displays plateaus in the localized phase; this feature 
is not only robust to the presence of a confining potential, but also to the presence 
of a finite entropy. In particular we observe that the overlap is most sensitive to the 
localization transition when k is close to the zone center, because this is the wave- 
vector region in which the features in S(k,oj) display the strongest k dependence in 
the delocalized phase. Unsurprisingly, the sharpest features in O at the localization 
transition are displayed for the largest Ak separation we investigated. Therefore the 
overlap provides a quantitative estimate of the localization transition based on the 
dynamic structure factor. Its effectiveness relies upon the fact that, in the problem at 
hand, all excitation modes localize and lose their dispersive nature at the ground-state 
transition point V = 4J. Yet it can be extended to situations (as the one discussed 
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in Section 4.1) in which only a part of the excited states in the excitation spectrum 
localizes. Most importantly, the overlap is a convenient experimental observable which 
allows to detect the transition even when one does not have access experimentally to 
the whole (k, u) plane. Indeed several recent experiments [6U [36j [39j SO] have measured 
w-scans for a fixed A; or a small range of k values (limited by the optical access to the 
sample). We observe that the overlap is well adapted to such a situation, providing a 
direct inspection into the localization transition by using only two tomographic scans 
of the dynamic structure factor at near wave-vectors. 



6. Summary and conclusions 

In conclusion, we have shown that the study of the dynamic structure factor provides 
important information on the interplay of disorder and interactions in one dimensional 
Bose fluids. Our analysis relies upon exact diagonalization results at arbitrary 
interactions, Bogolyubov theory for the weakly interacting case, and the exact solution 
for the hardcore case. Already in the clean case, the dynamic structure factor displays 
different features in the various interaction regimes, from sharp dispersive features in 
the weakly interacting superfluid phase, to a broad particle-hole continuum in the 
strongly-interacting (but still superfluid) phase, coexisting with a doublon mode. Once 
the disorder is turned on, we have shown that the dynamic structure factor allows to 
explore the spatial support of the excitations, yielding information on their localization 
properties. We have also investigated the features of the dynamic structure factor 
in the presence of a random-box disorder, from weak to strong interactions, showing 
that the dynamic structure factor captures the differences in the spectral features of the 
excitations with respect to the quasi-periodic potential. Thereby, it can probe in a direct 
way the localization mechanism. Finally, by exploiting the exact solution available for 
the hardcore-boson limit, we have shown that the main features of localization exhibited 
by the dynamic structure factor at zero temperature remain visible at a realistic finite 
entropy; and we have suggested an experimentally viable method to extract information 
on the Bose-glass transition at strong interactions, by analyzing frequency scans in the 
dynamic structure factor at two fixed and near wave vectors. While our paper exclusively 
focuses on the one- dimensional case, we argue that the insight provided by the dynamic 
structure factor into the physics of disordered bosons will be extremely useful also in 
higher- dimensional cases. In particular, the sensitivity of the dynamic structure factor 
to the extended or localized nature of the excitation modes, in a given frequency range, 
makes it a most viable probe of the presence of a mobility edge in the spectrum, which 
is a characteristic feature of higher-dimensional systems E| . 
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